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’9  ABSTRACT  {( 


Tomographic  visualization  of  three-dimensional  flow  fields  from  multidirectional 
interferometric  data  has  various  advantages  in  its  noninvasiveness  and  relatively  high 
resolution  and  instantaneous  capture  of  gross  fields.  The  technique,  however, 
encounters  ill-posed  problems:  that  is,  incomplete  projection,  limited  angular  scanning, 
and  insufficient  nonuniform  data.  In  the  research,  three  typical  methods  have  been 
developed  for  accurately  reconstructing  flow  fields  from  severely  limited  data.  First,  a 
general  method  termed  the  complementary  field  method  (CFM)  has  been  developed  in 
order  to  treat  the  ail  ill-posed  problems  in  a  unified  manner.  This  iterative  method  can 
be  combined  with  any  direct  reconstruction  techniques.  A  special  approach,  which 
employs  the  CFM,  has  been  formulated  for  discontinuous  shock  reconstruction  arising 
in  aerodynamics.  Second,  direct  reconstruction  techniques  based  on  continuous  local 
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v  basis  functions  have  been  developed.  These  techniques,  utilizing  higher  order 
approximation,  demonstrate  better  reconstruction  for  continuous  fields.  Third,  a  variable 
grid  technique  has  been  developed  to  reflect  intrinsic  spatial  resolution  information 
contained  in  interferometric  data.  Test  results  demonstrate  substantial  error  reduction 
when  the  CFM  and  an  appropriate  developed  technique  are  coupled,  f . 
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1. 


Statement  of  Problems  Investigated 


Interferometric  tomography,  that  is,  reconstruction  of  three-dimensional  refractive- 
index  fields  from  multidirectional  optical  pathlength  integrals,  has  great  potential  as  a 
powerful  diagnostic  tool  in  high-speed  aerodynamics.  The  interferometric  method  is  of 
great  value  in  its  instantaneous  capture  of  a  gross  field  and  noninvasiveness.  In  the  past, 
the  interferometric  reconstruction  has  been  implemented  by  those  methods  in  x-ray 
tomography;  however,  characteristics  of  these  two  fields  are  distinct. 


The  fields  in  high-speed  aerodynamics  are  continuous  except  at  a  shock.  Being 
sampled  usually  at  bright  and  dark  fringes,  the  original  interferometric  data  are  relatively 
sparse  and  nonuniform.  Especially,  near  a  shock,  the  sudden  density  change 
accompanies  data  loss  due  to  high  fringe  density  and  refraction.  In  addition, 
interferometric  reconstruction  frequently  encounter  major  difficulties  because  of  severely 
limited  projection  data  in  most  applications.  In  fluid  flow  measurements,  physical 
constraints  due  to  the  experimental  setup,  which  include  test  models  and  test-section 
enclosures,  are  frequently  present.  These  restrict  angular  scanning  or  block  a  portion  of 
the  probing  rays.  The  problems  of  incomplete  projection  and  limited  view  angle  thus 
occur.  The  adverse  effects  prohibit  accurate  reconstruction  of  object  fields  and  thus  limit 
practical  applications.  Accurate  numerical  reconstruction  requires  full  in'o.mation.  The 
ill-posed  data  reconstruction  is  sensitive  to  noise  and  produces  geometric  distortions  with 


various  artifacts.  In  order  to  accurately  reconstruct  flow  fields  from  ill-posed  conditions, 
new  techniques  need  to  be  developed,  which  incorporate  the  distinct  features  of  the 


interferometric  measurements. 


□  □ 


2.  Summary  of  Important  Results 

In  the  contract  period,  three  typical  new  techniques  were  developed  as  outlined 
below.  As  mentioned,  these  techniques  reflected  peculiar  aspects  of  interferometric 
tomographic  flow  measurements. 

2.1  Complementary  Field  Method 

A  general  reconstruction  method,  termed  the  complementary  field  method  (CFM), 
was  developed  to  deal  with  all  ill-posed  data  problems  in  a  unified  manner.  The 
developed  method  is  based  on  successive  estimation  of  the  flow  field  and  reconstruction 
of  the  corresponding  complementary  field  which  is  the  difference  between  the  estimate 
and  the  exact  field.  Since  the  complementary  field,  which  eventually  converges  to  a  null 
field  in  iteration,  is  reconstructed  instead  of  the  flow  field,  it  can  provide  better 
reconstruction  properties  and  less  error  energies.  Imperfect  inversion  techniques,  in 
addition  to  ill-posed  data,  contribute  to  reconstruction  errors.  The  CFM  can  be  employed 
to  abate  deficiencies  associated  with  both  measurement  data  and  compr  ational 
reconstruction  techniques.  Being  iterative  in  mature,  it  can  easily  incorporate  a  priori 
information.  In  reconstructing  a  flow  field,  sometimes  approximate  information  can  be 
known;  i.e.,  field  shape,  range  of  values,  discontinuity,  boundary,  etc,.  Even  slight 
knowledge  of  the  field  or  its  transform  can  be  efficiently  utilized  with  the  method.  The 
developed  method  was  tested  through  computer  simulation  of  experiments,  that  is, 
numerical  generation  of  ill-posed  data  and  corresponding  reconstruction.  When  no 
estimation  of  the  fields  is  made  in  the  initial  iteration,  the  fields  reconstructed  by  the  CFM 
displayed  gradual  reductions  in  their  maximum  as  well  as  average  errors,  which 
approached  limiting  values.  On  the  contrary,  the  conventional  methods  reached  an 
optimum  estimate  and  then  diverged.  The  test  observations  indicated  the  importance  of 
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the  initial  estimate  in  reducing  reconstruction  errors  and  expediting  convergence. 
Maximum  error  reduction  was  substantial.  Noticeable  improvements  occurred  near  the 
region  where  the  fields  obtained  a  local  maximum.  In  real  application,  this  region  is 
usually  of  our  interest.  The  test  results  of  simple  continuous  fields  were  encouraging. 
The  convergence  of  the  method  was  stable  but  rather  slow.  It  can  be  expedited  by 
employing  a  large  correction  parameter  in  calculating  a  new  estimate. 

2.2  Shock  Reconstruction 

Accurate  reconstruction  of  a  shock  requires  dense  sampling  at  the  discontinuity 
in  order  to  obtain  high  frequency  components  of  the  spatial  variation.  Due  to  the  data 
loss  or  inaccuracy  near  a  shock,  as  explained,  interferograms  may  not  produce  adequate 
information  to  resolve  the  shock  vicinity.  Consequently,  conventional  reconstruction 
techniques  are  not  appropriate  for  reconstruction  of  continuous  flow  fields  involving  a 
discontinuous  shock.  A  new  approach  was  formulated,  which  utilized  the  CFM  to 
improve  reconstruction  accuracy  of  a  locally  discontinuous  field.  The  method  is  based 
on  the  use  of  a  priori  information  of  the  shock  position  to  compensate  for  the  deficiencies 
associated  with  computational  techniques  in  reconstructing  a  discontinuous  field  and  loss 
of  information  due  to  noise.  In  a  test  of  a  high-speed  aerodynamic  flow,  the  shock 
positions  can  be  accurately  assessed  by  surface-mounted  transducers  or  other  optical 
means.  The  computer  simulation  of  experiments  indicated  that  fairly  accurate 
reconstruction  can  be  achieved  by  the  developed  method.  Error  reduction  was 
substantial  in  the  badly  defined  shock  vicinity  as  compared  with  conventional  methods. 
The  reconstruction  resulted  in  fast  stable  convergence. 
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2.3  Continuous  Local  Basis  Functions 


Conventional  methods  have  difficulties  in  accurately  reconstructing  a  continuous 
flow  field  from  sparsely  and  nonuniformly  sampled  interferometric  data.  Continuous  fields 
can  be  best  reconstructed  by  higher-order  continuous  local  basis  functions  (CLBFs). 
Series-expansion  techniques  using  higher-order  CLBFs  were  developed  for  reconstruction 
of  continuous  fields.  CLBFs,  having  support  in  a  local  region,  allow  precise  continuous 
approximation  of  a  field.  The  developed  techniques  demonstrated  power  of 
approximation,  independence  of  basis  functions,  and  efficiency  in  computing  projection 
data.  Power  of  approximation  provides  approximating  accuracy  and  smoothness  due  to 
the  continuity  of  basis  functions  up  to  higher-order  derivatives.  Overlapping  of  CLBFs 
ensures  a  sufficient  degree  of  freedom  to  represent  a  field  in  a  basis.  The  strong 
independence  of  the  overlapping  basis  functions  allows  uniqueness  in  determining  the 
coefficients  in  series  expansion.  It  guarantees  that  a  distinct  pair  of  coefficients 
corresponds  to  a  significantly  different  reconstruction.  Efficiency  for  projection  calculation 
enables  repeated  calculation  of  projection  integrals  of  CLBFs  in  setting  up  linear  algebraic 
equations  to  be  accomplished  with  great  simplicity.  The  CLBFs  allow  high  reconstruction 
accuracy  with  a  fewer  number  of  series  terms  and  minimal  artifacts  under  ill-posed 
conditions  and  undersampling.  Three  types  of  CLBFs  were  formulated  and  tested.  Being 
series  expansion  methods,  they  are  applicable,  directly  or  indirectly,  to  a  broad  range  of 
ill-posed  problems.  Linder  full  information,  performances  of  all  methods  (conventional 
and  CLBFs)  were  fairly  comparable.  However,  as  the  degree  of  ill-posedness  increased, 
the  CLBF  method  outperformed  conventional  ones,  especially  the  fixed  grid  method.  The 
fixed  grid  method  exhibits  large  errors  despite  the  post-smoothing  through  interpolation. 
When  noise  was  introduced,  reconstruction  errors  increased  with  noise  levels  in  all 
methods;  however,  the  CLBFs  were  less  susceptible  to  noise  than  other  methods.  The 
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CLBFs  could  be  defined  only  in  the  region  where  information  is  available.  They  could 
thus  produce  a  better-conditioned  coefficient  matrix  than  nonlocal  basis  functions, 
especially  under  the  presence  of  a  large  opaque  object. 

The  CLBF  method  was  combined  with  the  CFM  and  the  reconstructions  were 
compared  with  those  from  the  convolution  backprojection  and  CFM  combination.  Both 
combinations  reduced  reconstruction  errors  with  the  CFM  incorporated.  Under  full 
information,  reconstructions  were  comparable.  The  performance  of  the  CLBF 
combination,  however,  became  superior  as  the  degree  of  limited  data  and  the  level  of 
noise  increased.  Under  severely  limited-data  conditions,  the  CFM  reconstructions  initially 
converged  and  then  slowly  diverged. 

Variable  Grid  Method 

Interferometric  data  under  ill-posed  conditions  lead  to  nonuniform  scanning  of  a 
field.  This  nonuniform  scanning  should  provide  nonuniform  reconstruction  resolution  with 
greater  accuracy  in  the  region  that  rays  scan  more  closely.  In  general,  we  are  interested 
in  the  region  with  a  large  density  gradient.  This  region,  producing  dense  fringes,  allow 
close  sampling.  In  practice,  some  regions  can  be  heavily  scanned  to  extract  more 
information.  When  the  fixed  grid  method  (FGM,  a  conventional  method)is  employed,  a 
relatively  large  number  of  uniform  data  is  required  in  order  for  individual  grid  elements  to 
have  a  sufficient  number  of  rays  crossing  through  them.  In  this  way,  a  well-behaved 
coefficient  matrix  can  be  formulated  but  it  restricts  the  number  of  grid  elements, 
especially  under  ill-posed  conditions.  This  strategy  thus  limits  the  reconstruction 


resolution. 


In  the  research,  a  series  expansion  method  termed  variable  grid  method  (VGM) 
was  developed  to  reflect  the  intrinsic  spatial  resolution  information  contained  in  the 
original  measured  data  while  producing  a  well-behaved  coefficient  matrix  for  a  stable 
solution.  The  VGM  employs  nonuniform  rectangular-pulse  local  basis  functions.  In  the 
method,  the  bases  of  a  field  are  adjusted  so  as  an  appropriate  number  of  rays  pass 
through  each  of  them.  This  allows  a  reasonably  uniform  degree  of  redundancy  to  an 
individual  coefficient  in  the  linear  algebraic  equations.  The  resulting  well-behaved  matrix 
yields  a  stable  solution  of  the  coefficients.  The  algorithm  enables  the  region  of  interest, 
scanned  densely  with  probing  rays,  to  contain  fine  local  bases  and  thus  to  be 
reconstructed  with  higher  resolution. 

The  following  general  interpretations  can  be  drawn  on  the  reconstruction  result 
comparison  between  the  VGM  and  FGM.  For  complete  data,  the  performance  of  the  two 
methods,  including  the  number  of  bases,  was  fairly  comparable.  They  produced  the 
almost  same  number  of  the  optimum  bases.  As  the  degree  of  ill  cosedness  (incomplete 
projection  and  limited  view  angle)  increased,  especially  limited  view  angle,  the  VGM 
produced  a  large  number  of  appropriately  adjusted  bases.  This  led  finer  resolution  in  the 
region  of  interest  and  thus  less  reconstruction  errors.  Maximum  errors  were  reduced 
mostly  in  the  areas  with  large  field  values  or  steep  gradients.  On  the  contrary,  the  FGM 
produced  appreciable  distortions  in  reconstruction.  This  is  believed  due  to  the  ill-behaved 
coefficient  matrix.  Overall  the  VGM  performance  was  better  than  that  of  the  FGM  under 
any  circumstances.  The  VGM  allowed  higher  reconstruction  resolution  and  accuracy  for 
a  given  number  of  series  terms.  It  thus  produced  minimal  artifacts  even  under  ill-posed 
conditions  and  undersampling. 
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Finally,  the  VGM  can  be  combined  with  the  CFM  for  further  performance 
enhancement.  This  combination  needs  to  be  tested  and  compared  with  other 
combinations. 
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